Orientations, lattice polytopes, and group arrangements II: 
Modular and integral flow polynomials of graphs 

Beifang Chen and Richard Stanley 

Abstract. We study modular and integral flow polynomials of graphs by 
means of subgroup arrangements and lattice polytopes. We introduce an Euler- 
ian equivalence relation on orientations, flow arrangements, and flow poly- 
topes; and we apply the theory of Ehrhart polynomials to obtain properties 
of modular and integral flow polynomials. The emphasis is on the geometri- 
cal treatment through subgroup arrangements and Ehrhart polynomials. Such 
viewpoint leads to a reciprocity law for the modular flow polynomial, which 
gives rise to an interpretation on the values of the modular flow polynomial at 
negative integers, and answers a question by Beck and Zaslavsky. 



1. Introduction 

The flow polynomial ^{G, t) of a graph G was introduced by Tutte [24] as a 
conceptual dual to the chromatic (or tension) polynomial of G. When G is a planar 
graph, tf{G, t) is essentially the chromatic polynomial i) of the dual graph G* 

in the sense that x{G* ,t) = t'^'^'^^ ip{G,t), where c{G) is the number of connected 
components of G. The historic Four-Color Conjecture of the time was made by 
Tuttc into the Five-Flow Conjecture: any bridgeless graph admits a nowhere-zero 
integer 5-flow. Both conjectures are still open and stimulate studies on chromatic 
and how polynomials. In a seminal paper |19| . Rota introduced characteristic 
polynomial for posets and observed that ^p{G, t) is the characteristic polynomial 
of the circuit lattice of G. Greene and Zaslavsky |15) made Rota's observation 
transparent between the flow polynomial and the characteristic polynomial by using 
hyperplane arrangements. As a special case of Zaslavsky's formula [25 , the absolute 
value \'p{G, — 1)| counts the number of totally cyclic orientations of G, which is a 
dual analog of Stanley's result on chromatic polynomials: |x(G, — 1)| counts the 
number of acyclic orientations of G. However, Stanley's result |20j includes an 
interpretation of the values of x(G, t) at negative integers, known as the Reciprocity 
Law of chromatic polynomials. More recently, Kochol [17] showed that the number 
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of integer- valued g-flows is a polynomial function of q and introduced the integral 
flow polynomial (p^{G,t); Beck and Zaslavsky |1] studied the modular and integral 
flow polynomials for graphs and signed graphs, using Ehrhart polynomials of lattice 
polytopes. The present paper, as a continuation of [11 , is to associate flow group 
arrangements with graphs, and to obtain a clear picture of the relation between 
the integral flow polynomial and the modular flow polynomial. The byproduct of 
this association is a generalization of the Reciprocity Law of chromatic and tension 
polynomials to modular and integral flow polynomials, and the interpretation of the 
values of the modular and integral flow polynomials at zero and negative integers. 
The geometric method of our exposition may be modified to obtain analogous 
results on Tutte polynomials. 

Let G = (y, E) be a finite graph with possible loops and multiple edges. We 
write V = V{G), E = E{G). For each subset X <Z E, denote by {X) the induced 
subgraph {V,X). An orientation on G is a (multivalued) function e : V x E ^ 
{—1, 0, 1} such that (i) e{v, e) has the ordered double- value ±1 or =f1 if e is a loop at 
a vertex v and has a single- value otherwise, (ii) e(u, e) = if w is not an end-vertex 
of e, and (iii) e{u, e)e{v, e) = —1 if e has two distinct end-vertices u,v. Pictorially, 
if e is a non-loop edge with distinct end-vertices u,v, then e{u,e) = —e{v,e) = 1 
(or e(v, e) = —£{u, e) = 1), and it means that e is assigned an arrow from u to v, 
which contributes exactly one out-degree at u and one in-degree at v. If e is a loop 
at a vertex v, then e{v, e) = ±1 or =pl, and it means that the loop e is assigned an 
arrow, pointing away and to the vertex v, which contributes exactly one out-degree 
and one in-degree at v. We assume —(±1) = T^, ^ (t1) = il- A graph G together 
with an orientation e is called a digraph^ denoted (G, e). A digraph is said to be 
directed Eulerian if its in-degree equals its out-degree at every vertex. 

Let (G, e) be a digraph throughout the whole paper. Associated with (G, e) is 
the incidence matrix M = M{G) :— [niy^ejvxE, where rriy^e = if the edge e is 
a loop and m^j g = e{v, e) if e is not a loop. Let A be an abelian group. A flow 
of (G, e) with values in A, or an A-flow, is a function f : E A, satisfying the 
Conservation Law: 

^m,,e(«,e)/(e) =0 or ^ e(z;, e)/(e) = 0, v e V, (l.I) 

eeE e£E 

where e(w, e) is counted twice in the second sum as —1 and 1 if e is a loop at 
its unique end-vertex v. A flow / is said to be nowhere-zero if /(e) ^ for all 
eeE. We denote by F{G,e]A) the abelian group of aU A-flows of {G,e). The 
flow arrangement of (G, e) with the abelian group A is the group arrangement 
Afl{G, e; A) of F{G, e; A), consisting of the subgroups 

F,:={f eF{G,e;A)\f{e)=0}, eeE. (1.2) 

Coincidentally, we shall see that the characteristic polynomial x{-^Fh{G, e; A),t) is 
equal to the (modular) flow polynomial (p{G,t), defined for t = q as the number 
of nowhere-zero flows of (G, e) with values in an abelian group of order q. The 
polynomial ip is independent of the chosen orientation e and the abelian group 
structure; see Rota |19) and Tutte |23) . For a complete information about modular 
and integral flows, we refer to the book of Zhang |27| . 

Recall that a cut of G is a nonempty edge subset of the form [S, S'^] , where 
S <e V is a nonempty proper subset, S'^ := V — S is the complement of S, and 
[S, S'^] is the set of all edges between the vertices of S and the vertices of 5"^. Let 
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U — [S, S'^] be a cut. A direction of U is an orientation su on the induced subgraph 
(y, U) such that the arrows of the edges in U are either all from S to S'^ or all from 
S'^ to S; U together with a direction eu is called a directed cut, denoted {U,eu)- 
If {U, e) is a directed cut, we call {U, e) a directed cut of both (G, e) and e, and 
say that the cut U is directed in {G,e). Let 0(0) denote the set of all orientations 
on G. We denote by Otc{G) the set of all orientations without directed cut (also 
known as totally cyclic orientations, as they are the orientations in which every 
edge belongs to a directed circuit). 

Let (pz{G, q) denote the number of nowhere-zero integer- valued flows / of (G, e) 
such that < |/(e)| < q for all e e -E. As pointed out by Beck and Zaslavsky pQ, 
the function ip^ (G, q) was never mentioned to be a polynomial until Kochol [17j . 
If e is totally cyclic, we introduce the counting functions 

MG,q) := #{./ e F{G,e;Z) \ < /(e) <q,ee E}, (1.3) 

Ve{G,q) := #{/ e F{G,e;Z) \ < /(e) < g, e G E}, (1.4) 

and the relatively open 0-1 polytope 

A+ (G, e) := {/ G F(G, e; R) | < /(e) < 1, e G E}. (1.5) 

The closure A+(G, e) is a 0-1 polytope (whose vertices are 0-1 vectors), and is the 
convex hull of all 0-1 flows of (G, e) in M.^ (flows whose values are either or 1), 
called the flow polytope of (G, e). We introduce the following counting function 

^.(G, q) : = #{(p, /) I p G Otc(G), / G i^(G, p; Z), < /(e) < g, e G E}. (1.6) 

We shall see that ^Pz{G, q), fe{G, q) are polynomial functions of positive integers q, 
and (pz{G, q), (pe{G, q) are polynomial functions of nonnegative integers q, and that 
(p^{G, q) is independent of the chosen orientation e. The corresponding polynomial 
ip^{G,t) {(p^{G,t)) is called the (dual) integral flow polynomial of G, and tpe{G,t) 
{(pe{G,t)) the local (dual) flow polynomial with respect to the orientation e. The 
names and notations are so selected in order to easily recognize these polynomials. 

We first reproduce a result due to Kochol '17 about Equation (|1.8|) . and due 
to Beck and Zaslavsky 1 about the combinatorial interpretation of the values of 
ip^{G,t) at nonpositive integers. 

Theorem 1.1 (Kochol [17 , Beck and Zaslavsky [1 ). Let G = {V,E) he a finite 
bridgeless graph with possible loops and multiple edges. 

(a) // the orientation e is totally cyclic, then /S.'^^{G,£) is a relatively open 
0-1 polytope in M.^ of dimension n[G); Lpe{G,t) and (pe{G,t) are Ehrhart 
polynomials o/ A+(G, e) and A+ (G,£) respectively, and satisfy the Reci- 
procity Law: 

^,(G,-0 = (-l)"('=^^e(G,t), (1.7) 

where n{G) = \E\ — r{G) and r{G) is the number of edges of a maximal 
spanning forest of G. Moreover, 



^e(G,0) = (-l)"(«), ^,(G,0) = 1. 
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(b) The integral flow polynomials ip^{G,t) and ip^{G,t) can be written as 

MG,t) = VpiG.t), (1.8) 

peOTc(G) 

<^.(G,t) = ^piG,t). (1.9) 

pGOTr(G) 

and satisfy the Reciprocity Law: 

(/..(G,-0-(-l)"(^^^.(G,t). (1.10) 

In particular, \ip^{G,0)\ counts the number of totally cyclic orientations 
on G. 

There are analogous results on the modular flow polynomial Lp{G,t). To do 
this we need to introduce an equivalence relation on the set 0{G) of orientations 
on G. Two orientations £1,62 on G are said to be Eulerian equivalent, written 
£1 ^ £2, if the spanning subgraph induced by the edge subset {e ^ E \ £i(u, e) ^ 
£2(w,e)} is a directed Eulerian subgraph with respect to the orientation either £1 
or £2. We shall see that ~ is indeed an equivalence relation on 0{G). Moreover, 
if an Eulerian equivalence class intersects C'tc(G), the whole equivalence class is 
contained in C'tc(G). So ^ induces an equivalence relation on the set C'tc(G) of 
totally cyclic orientations. Let [(!?tc(G)] denote a set of distinct representatives, 
exact one representative from each equivalence class of ~ on Otc(G). We introduce 
the following counting function 

^(G, q) : = #{(p, /) I p e [Otc(G)], / G F{G, p; Z), < /(e) < g, e e E]. (1.11) 

We next produce the following Theorem 1 1.2[ which answers a question by Beck and 
Zaslavsky [1] about the combinatorial interpretation of the values of the modular 
flow polynomial ip{G,t) at zero and negative integers. 

Theorem 1.2. Let G = {V,E) be a finite bridgeless graph with possible loops and 
multiple edges. Then ip(G,q) (ip{G,q)) is a polynomial function of degree n{G) of 
positive (nonnegative) integers q, and satisfy the Reciprocity Law: 

^(G,-t) = (-l)«(«)^(G,t). (1.12) 

Moreover, 

if{G,t)= MG,t), (1.13) 

pe[0-rc(G)] 

^{G,t)^ Y ^p{G,t). (1.14) 

pe[OTc(G)] 

In particular, \ip{G,—\)\ counts the number of totally cyclic orientations on G, 
and \ip{G,Q)\ counts the number of Eulerian equivalence classes of totally cyclic 
orientations. 

Equation (|1.13p is recently obtained by Kochol ^17, with a formal proof in 
different form. The combinatorial interpretation of |(/3(G, — 1)| = rG'(0, 2) is due 
to Las Vergnas |18) . see also Brylawski and Oxley [6]. At the moment of this 
revising, we noticed a paper by Breuer and Sanyal [5] on modular flow reciprocity, 
which is quite different from our Reciprocity Law (|1.12p . The difference lies in that 
the result of [5] on i^(G, —q) for a positive integer q involves the counting of flows 
modulo q, our result only involves nonnegative integer flows bounded by g, and 
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that the bijection between the two counting sets is nontrivial; see Section 5 for the 
detailed discussion. The proof of Theorem 11.21 is rigorous and self-contained. The 
following corollary is an immediate consequence of Theorem 11.21 

Corollary 1.3. The value Tg(0,1) of the Tutte polynomial TG{x,y) counts the 
number of Eulerian equivalence classes of totally cyclic orientations on G. 

2. Characteristic polynomials of group arrangements 

Let f2 be a finitely generated abelian group. By a flat of we mean a coset 
of a subgroup of f2. For a subgroup F C f], we denote by Tor(F) the torsion sub- 
group of F and write |F| := |Tor(F)| ^■■ank(r)^ gy ^ subgroup arrangement (or just 
arrangement) of we mean a finite collection of flats of fi. Associated with a sub- 
group arrangement A is the semilattice ^{A), whose members are nonempty sets 
obtained from all possible intersections of flats in A. The characteristic polynomial 
of A is defined as 

where /i is the Mobius function of the poset ^{A), whose partial order is the set 
inclusion, {X) := {x — y \ x,y ^ X} . 

Let ^(ri) be the Boolean algebra generated by cosets of all subgroups of 51, 
i.e., every member of ^(Sl) is obtained from cosets of subgroups of VI by taking 
unions, intersections, and complements finitely many times. A valuation on with 
values in an abelian group A is a, map v : £s§{Q) — >■ A such that 

1^(0) -0, 

i^{X UY) = i^{X) + i^{Y) - v{X n Y) 
for X,Y ^ ^(fl). A valuation is said to be translation invariant if 

iy{S + x) = iy{S) 

for S G ■S§{^) and any a; e fi; and v is said to satisfy multiplicativity if 

v{A + B) = iy{A)iy{B) 

for subgroups A, B C such that A -f _B is a direct sum of A and B, and the 
subgroup j4 + _B is a direct summand of fl. 

Theorem 2.1 (Chen [11] ). For any finitely generated abelian group f2, there exists 
a unique translation invariant valuation A : £s§{Q) — > Q[t] such that the multiplica- 
tivity is satisfied and 

x{n) = |Tor(i7)|i™^'^(^) = 

In particular, X{T) = Jnjri '^^V subgroup F C and for any subgroup arrange- 
ment A of Q, 

x(n~ y x)^xiA,t). 

The analogue of Theorem l2.1l for vector spaces was obtained by Ehrenborg and 
Readdy |13j . Let y be a vector space over an infinite field. Let ^(V) be the lattice 
of all affine subspaces of V . We denote by SSiy ) the Boolean algebra generated by 
I£{y ). A subspace arrangement of V is a finite collection A of affine subspaces of V. 
The torsion of any subspace is just the zero space. Then characteristic polynomial 
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t) of a subspace arrangement A can be defined by the same formula p.ip for 
subgroup arrangement. 

Theorem 2.2 (Ehrenborg and Readdy 13 ). For any finite- dimensional vector 
space V over an infinite field IK, there exists a unique translation invariant valuation 
A : ^{V) -J- Z[t] such that X{W) = i'^™^ for subspaces W CV. Moreover, for a 
subspace arrangement A of V , 

^(^- U ^) -=x(Ai)- 

^ XeA ' 

One may combine Theorems 12.11 and 12.21 by considering arrangements of afhne 
submodules. Let M be a finitely generated left i?-module over a commutative ring 
i?; we restrict R to the cases of M, Z, and Ij/qL. By a flat of M we mean a subset 
of the form a-\-N = {a-\-x\x^ A^}, where is a submodule of M . Let ^{M) be 
the lattice of all flats of M, and 3§{M) the Boolean algebra generated by ^{M). 
For each subset S C M, we denote by I5 the characteristic function of S. 

Let „4 be a finite collection of flats in M , called a submodule arrangement of M. 
Let jSf (.4) be the poset whose members are nonempty sets obtained by taking all 
possible intersections of flats in A, and whose partial order < is the set inclusion. 
For each X e J^f{A), we define 

X° —X - [j Y. 

Ye^(A), Y<X 

Clearly, {X° \ X e ^{A)} is a family of disjoint subsets of M. Then for each 

Xe^{A), 

Ix = ^ 

Ye^(A),Y<X 

By the Mobius inversion, for each X e ^{A), 

Ye^(A), Y<X 

In particular, M° = M -{jA = M - [jxeA ^ 

1a/-U^= E l^i^^M)lY. (2.2) 

Yese(A) 

Thus for any valuation v on B{M), we have the Inclusion- Exclusion Formula: 

iy{M-[jA)= KX,M)v{X). (2.3) 

xeJ^iA) 

This is a prototype of many existing formulas when v is taken to be various valua- 
tions; see P [131126] . 

Let M be the Euclidean n-space R". One has half-spaces {x G R" \ f{x) < c} 
(with linear functionals / : R" R and constant real numbers c) , convex polyhedra 
(intersections of half-spaces), and the Boolean algebra ^("P") (generated by half- 
spaces by taking intersections, unions, and relatively complements finitely many 
times). There are two valuations x ^-nd x on 3§{'P"), both are referred to the 
Euler characteristic (see O 1221 I26j . for example), such that for any relatively 
open convex polyhedron P, 

X{P) = (-l)'^™^, x{P) = lim x(^n [-r,r]"). 
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By Groemer's extension theorem |16j . x x can be extended to be hnear func- 
tionals on the functional space spanned by characteristic functions of convex poly- 
hedra. Now let ^ be a hyperplane arrangement of M". Evaluating x ^^id x on 
both sides of (|2.3p , one obtains Zaslavsky's first and second counting formulas (see 
pl[25l[26] ): 

|x(^, -1)1 = number of regions of R" - IJ-^: (2-4) 
|x(^, 1)1 = number of relatively bounded regions of R" — IJ^. (2-5) 

3. Modular flow polynomials 

Let {Hi, Si) be subdigraphs of the graph G = {V,E), i = 1,2. The coupling of 
El and £2 is a function [£1,62] : E — >• { — 1,0, 1}, defined for each edge e G (at its 
one end- vertex f) by 

r 1 if eeS(i/i)nS(i72),£i(",e)= £2(1^,6), 
[ei,e2\{e)^l -1 ii e £ E{H^) f^ EiH^), e^{v,e) ^ e2{v,e), (3.1) 
[ otherwise. 

The following proposition is straightforward. 

Proposition 3.1 (Berge [2]). (a) ^ function f : E{G) A is a flow of 

(G,e) if and only if for any directed cut {U,eu), 

J2heu]{e)f{e)^0. 

eeu 

(b) In particular, the digraph {G,e) is directed Eulerian if and only if for any 
directed cut {U,eu), 

^[e,ec/](e) = 0. 

(c) The graph G is Eulerian if and only if every cut U contains even number 
of edges. 

Let Fnz(G,£; A) denote the set of all nowhere-zero flows with values in A, i.e., 

Fn,(G, e- A) := {/ e F(G, e; A) \ /(e) ^ 0, e e E}. 

If |yl| — q\s finite, it is well-known (see [27]) that the counting function 

^(G,g) := |i^„,(G,£;A)| (3.2) 

is a polynomial function of g, depending only on the order \A\, but not on the 
chosen orientation e and the group structure of A. The polynomial <p(G, t) is called 
the modular flow polynomial of G. 

For two orientations p,cr ^ C'(G), there is an involution Pp^a '■ A^ A^ , 
defined by 

^ ^ \ -/(e) li p{v,e) ^ a{v,e). ^ ' 

In fact, Pp^ef — [PiE]/- Obviously, Pp^p is the identity map, Pp^aP<j,e ~ Pp,e- 
Lemma 3.2. The involution Pp^^ is a group isomorphism. Moreover, 

Pp,,FiG,e;A) = F{G,p-A), 
Pp,eFUG,e;A) = F,,,{G,p-A). 
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Proof. It is clear that Pp^e is invertible and Ppl = Pp.e- Let / e . The 
group isomorphism follows from the fact that at each vertex w, 

^ p{v, e)(Pp,e/)(e) = ^ e)p{v, e)e{v, e)/(e) = ^ e{v, e)/(e). 

eeB eG-E eGB 

Since Ppef{e) 7^ is equivalent to /(e) ^ 0, it follows that Pp ^Fa.^{G,e; A) = 
F^,iG,p-A)- ' □ 

Let T be a maximal forest of G in the sense that every component of T is a 
spanning tree of a component of G. For each edge e of the complement T"^ := 
E — E{T), let Ce denote the unique circuit contained in T U e, and let pe be a 
direction of Ce (i-e. {Ce,Pe) is directed Eulerian) such that Pe(e) = £(e). It is easy 
to see that [e, pe] is a flow of (G, e). 

Lemma 3.3 (Berge 2 ). Let T be a maximal spanning forest of G. Then each flow 
f of the digraph (G, e) can he expressed as a unique linear combination 

f=Y.f{e)[e,pe]. (3.4) 

The system (jl.ip for a flow /, whose equations are indexed by vertices v V, 
is equivalent to the system p.4|) . which can be written as 

E /(e)[e,Pe](.T), xeT, (3.5) 

eGT<^ 

whose equations are indexed by edges x € T. In other words, the values /(e) for 
e £ T"^ can be arbitrarily specified, and f{x) for x G T are determined by p.Sp . 

Let n{G) denote the cycle rank of G. If T is a maximal spanning forest of G, 
then n(G) is the number of edges of T"^. Note that 

n(G) = 1^1 - IV^I +c(G), 

where c(G) is the number of connected components of G. Lemma 13.31 shows that 
the abelian group F{G,e;A) is of rank n{G). The flow arrangement of the di- 
graph (G, e) with the abelian group A is the subgroup arrangement Afl{G, e; A) of 
F{G, £; A), consisting of the subgroups 

Fe=FeiG,e;A):={f eF{G,e;A)\f{e)=0}, e e E. (3.6) 

The semilattice ^(^fl(G, e; A)) consists of the subgroups 

Fx = Fx{G, e; A) := {/ e F(G, e; A) \ f{x) = 0, x 6 X}, X C E. 

Then F = F0 = F{G,e;A). Using (13. 5|) . it is easy to see that the arrangement 
AphiG, e; A) is isomorphic to the arrangement Afl{G, T, e; A) of A^^'^''\ consisting 
of the subgroups 

{/ e A^(^') I /(e) =0} if e G 

{/ e A^(^') I E.GTc[£>P=^](e)/(^) - 0} if e 
The isomorphism is given by the restriction / f\E{T')^ sending Fe to He- 

Lemma 3.4. The abelian group F{G, e; A) is isomorphic to the product group 
^"(G) Moreover, for any subset X C E, the subgroup Fx{G,e;A) is isomorphic 
to the product group A^^-^^-^K In particular, if \A\ = q is finite, then 

\F{G,e;A)\^q''^^\ |Fx(G, e; A)| = (3.7) 



He 
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Proof. Trivial □ 

Theorem 3.5. Let A be an abelian group such that either \A\ — q is finite, or 
A = Tj, or a is an infinite field. Then 

^{G,q)^x(F{G,e;A)- \J 

^ eeE 

Moreover, \(p{G, — 1)| counts the number of totally cyclic orientations on G. 



^x{A,dG,e;A),q). (3.8) 

t=q 



Proof. Let X,Y C E. If \A\ > 2 (including \A\ = oo), then Fx{G,e;A) C 
Fy (G, e;A) is equivalent to that circuits of {E — X) are contained in {E~Y). Thus 
the map Fx{G,e; A) H> Fx{G,e;R) is an isomorphism from Jf{AFh{G,e; A)) to 
^(^fl(G, e; R)). Consequently, we have the same Mobius function 

tx{Fx (G, e; A), Fy (G, e;A))= fi{Fx (G, e; K), (G, e; R)). 

If A is infinite, applying the valuation A to both sides of (|2.2p . we have 



If A is finite and \A\ = q, applying the counting measure # to both sides of (|2.2p . 
we have 

The identity p.Sp follows immediately for positive integers q. 

For g = 1, we have A = {0} and F„z(G,e; A) = 0. Hence ip{G, 1) = 0. Since the 
hyperplane arrangement A^hiG, s;R) is central (all hyperplanes pass through the 
origin), there is no relatively bounded region. Zaslavsky's second counting formula 
((23)) confirms that x(Al(G,£;R), 1) = 0. 

Finally, by Zaslavsky's counting formula (j2.4p . (— l)"'^'')(y9(G, — 1) counts the 
number of regions of the complement F — Uee_E-^e- Lemmas 14.21 and 14.31 the 
regions of the complement F{G, e; M) — IJ Afl{G, e; M) correspond bijectively to the 
totally cyclic orientations on G. □ 

4. Integral flow polynomials 

In this section we apply the Ehrhart polynomial theory to study integral flow 
polynomials. Let us recall some well-known facts about lattice polytopes and 
Ehrhart polynomials. Let P be a relatively open lattice polytope of R", i.e., P 
is open in the flat that P spans, and the vertices of P are lattice points of Z". 
The closure of P is denoted by P. A bounded lattice polyhedron is a disjoint union 
of finitely many relatively open lattice polytopes. Let X be a bounded lattice 
polyhedron and q a positive integer. The dilatation of X by g is the polyhedron 
qX :— {qx \ x G X}. Let 

L{X,q) :=#(qXnZ"). 
It is known that L{X,q) is a polynomial function of degree dimX in the posi- 
tive integer variable g, called the Ehrhart polynomial of X. Moreover, the leading 
coefficient of L{X,t) is the volume of X; the constant term L{X,0) is the Euler 
characteristic xi-^)- In particular, if X = P is a relatively open lattice polytope, 
then L(P,q) and L[P,q) satisfy the Reciprocity Law: 

L{P,-t) = i-lf '^PL{P,t); 
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the constant term of L{P,t) is 1, and the constant term of L{P,t) is (— 1)"^'™^. 
All these and other related properties about Ehrhart polynomials can be found in 

[SIM mi- 

Flows with values in R are called real flows; and flows with values in Z are 
called integer flows. A flow / G F{G,e;M.) is called a q-flow if |/(e)| < q for all 
e € E. We define the set of all real g-flows of (G, e) as 

F{G,e;q) {/ £ F(G,e;M) : |/(e)| <q,ee E}. 

We denote by F^{G, e; q) the set of aU integer q-fiows of (G, e), and by Fnzz{G, e; q) 
the set of all nowhere-zero integer q-fiows, i.e., 

Fnz.(G, e; q) := {/ £ (G, e; g) | /(e) ^0,eeE}. 

Clearly, (G, e;q) is the set of lattice points of the dilatation gApL (G, e) (dilated 
by g) of the nonconvex polyhedron 

Ap,(G,£) := {/ e ^^(G,£;R) : < |/(e)| < l,e G -B}. 

It follows that the counting function 

MG,q) |Fnz.(G,£;g)| =L(Ap,(G,£),g) (4.1) 

is an Ehrhart polynomial of degree dim Afl(G, e) in the positive integer variable q. 
In fact, we shall see that \Fnzz{G, e; q)\ is independent of the chosen orientation e. 
We call ip^{G,t) the integral flow polynomial of G. 

Lemma 4.1. T/ie involution Pp^g is a group isomorphism from M.^ to itself. More- 
over, 

Pp,,Ap,(G,£) = Ap,(G,p), 
Pp,eFnz^{G,£;q) = Fnz^iG,p;q). 

Proof. Let A M, / e M^, ee E. Note that < |Pp,£/(e)| < 1 is equivalent 
to < |/(e)| < 1. Hence Pp^,A^^iG,e) = Apl(G,p). Similarly, Pp,,/(e) £ Z is 
equivalent to /(e) G Z; and < |Pp^e/(e)| < g is equivalent to < |/(e)| < q. Thus 
Pp,eFn.AG,e;q)^FnUG,p;q)- ' □ 

Let A = R. The subgroup arrangement .4fl(G, £;K) is a hyperplane arrange- 
ment of F{G, e; R). The complement of ^fl(G, £; R) is the set 

F{G, £; R) - U Al(G, £; R) = {/ G F{G, e; R) | /(e) ^ 0, e G E}. 

For each edge e G E with end-vertices m, w, the nonzero condition /(e) ^ can be 
split into two inequalities: 

/(e) >0 and /(e) < 0; 

the former can be interpreted as an orientation of e agreeing with s{u, e), and the 
latter is interpreted as an orientation of e opposite to e{u, e). 

For each orientation p G 0{G), we introduce the open convex cone 

GP{G,e) := {/ G P(G,£;R) : [p,e]{e)f{e) > 0,e G E}. 

The complement F{G, £;R) — IJ Al(G, e; R) is a disjoint union of these open convex 
cones, some of them may be empty. By Lemma [4.21 below, the cone C^{G,e) is 
isomorphic to the open convex cone 

G+(G,p) {/ G F(G,p;R) 1 /(e) > 0,e G E]. 
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We introduce the relatively open polytopes 

A+ (G, p) {/ e F{G, p; K) I < /(e) < 1, e e E}, 
A,^,(G,e) {/ e F(G,£;R) | < [p,e](e)/(e) < l,e G i?}. 
If A+ (G, /o) 7^ (equivalent to that p is totally cyclic) , then the closure of A+ (G, p) 
is the polytope 

A+ (G, p) {/ e ^^(G, p; R) I < /(e) < 1, e e E}. (4.2) 

Whether the orientation p is totally cyclic or not, the set A+(G, p) is always a 
polytope, and is called a /Zoiw polytope of G with respect to p. 

Lemma 4.2. Pp.jAFL(G,e) = A+(G,p); and disjoint decomposition 

Ap,(G,£) = □ A^,(G,£). (4.3) 

peo{G) 

Proof. Since Pp^ef = [p,e]/ for / G M^, then the first identity is trivial by 
definition of ApL(G,e), A+(G,/9), and Lemma 

Let / G Afl(G,£). We define an orientation p on G as follows: for each edge 
e at its one end-vertex w, set p{v,e) = e{v,e) if /(e) > and p{v,e) — —e{v,e) if 
/(e) < 0. Then / G Afl(G,£). Conversely, each Afl(G,£) is obviously contained 
in Afl(G, e). The union is clearly disjoint. □ 

Notice that the open convex cone G'^(G,e) may be empty for the given ori- 
entation e. If (G, e) contains a directed cut, then it is impossible to have positive 
real flows by Proposition OJa), thus C+(G,e) = 0. To have G+(G,e) ^ 0, the 
orientation s must be totally cyclic. 

Lemma 4.3. (a) A+ (G,e) 7^ «/ and only if e e C'tc(G). 

(b) If s & Otc(G), then A+(G,£) is a 0-1 polytope, i.e., all its vertices are 
flows of (G, e) with values in {0, 1}. 

(c) If s d 0Tc(G), then every flow of (G, e) with values in {0, 1} is a vertex of 
A+(G,£). 

Proof, (a) Trivial. 

(b) Let / be a vertex of A+ (G, e). It is enough to show that / is integer- valued. 
By Linear Programming the vertex (xe = /(e) : e G £') is a unique solution of a 
linear system of the form 

J2e£E'^v,(^^e = 0, veV, 
Xe = ae, e G Ef, 

where — or 1, Ej is an edge subset, \Ef \ = \E\ —n{G). The system is equivalent 
to the linear system 

where E'j^ :^ E — Ef, by G Z, and the rank of the matrix [mt,^e]yx_Ej. is n{G). 
Since the incidence matrix M — [rriy^elvxE is totally unimodular (see [4], p. 35), 
the submatrix [my elvxE'j, is row equivalent to the matrix [^j over Z, where / is 
the identity matrix when V is linearly labeled. It then follows that the solution 
{xe = /(e) : e G i?) is an integer vector. 

(c) Notice that a flow of (G,e) with values in {0, 1} is just the characteristic 
function of the edge set of a directed Eulerian subgraph of (G, e). Let / be a 
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flow with values in {0,1}. Suppose / is not a vertex of A+(G,£). Then there 



are distinct flows fi of {G,e) with values in {0, 1} such that / = X^iLi Oj/i, where 



> 0, X]i=i fli = 1, fc > 2. Let e be an edge such that /(e) = 1. Then fi{e) = 1 
for aU i; otherwise, say, /i(e) 0, then /(e) = ^ifM ^ J2i=2 < ^' which 

contradicts /(e) = 1. Thus fi — f for all i; this is contradictory to the distinctness 
of/.. □ 

Recall that Lp^{G,q) for a positive integer q is the number of integer flows of 
(G, e) with values in {1, 2, . . . , g — 1}. In other words, (pe{G, q) is the number of 
integer flows of (G,£) with values in the open interval (0,(j). Clearly, ip^{G,q) 
counts the number of lattice points of the dilatation gA+ (G,£), i.e., 

^aG,<z)=i(A+(G,£),<z). (4.4) 

We call ip^{G,q) the local flow polynomial of G with respect to e. Analogously, let 
(pe{G, q) denote the number of integer flows of (G, e) with values in {0, 1, . . . , g}. 
In other words, (pe{G,q) is the number of integer flows of (G, e) with values in 
the closed interval [0, g]. Then (ps{G,q) counts the number of lattice points of 
gA+(G,e), i.e., 

^aG,g) = L(A+(G,£),q). (4.5) 

We call (pe{G,q) the local dual flow polynomial of G with respect to e. Now we 
denote by (p^{G, q) the number of pairs (p, /), where p is a totally cyclic orientation 
on G and / is an integer flow of (G, p) with values in {0, 1, . . . , g}. We call <^zj(G, q) 
the dual integral flow polynomial of G. 

Proof of Theorem 11.11 

(a) By Lemma l473l the closure A+(G, e) of the open polytope A+(G, e) is the 
convex hull of the lattice points / £ such that /(e) G {0, 1} for all e G -E and 
satisfying (|l.ll) . The Reciprocity Law and the interpretation of the constant term 
follow from the Reciprocity Law and the properties of Ehrhart polynomials. 

(b) Note that Fnz^(G, £;(/) = g'AFL(G, e). By Lemma H^l we have a disjoint 
union 

(zAp,(G,£)= □ q/^PdG,e), (4.6) 

where each lattice open polytope ApL(G,e) is isomorphic to the 0-1 open polytope 
A+(G,p) by the unimodular transformation Pp^e- Then (11.81) follows immediately 
from (14. 6p : (|1.9p follows from deflnition of (p^{G,q). The Reciprocity Law (jl.lOp 
follows from (|1.7p - (ll.9l) . The interpretation of the constant term (p^{G,0) follows 
from dllHl) and (^e(G,0) (-1)"^'^'. □ 

5. Interpretation of modular flow polynomial 

This section is devoted to interpreting the values of the modular flow polynomial 
in a way similar to how the modular tension polynomial was interpreted in [llj . 
For the graph G ~ (V, E) and a positive integer q, there is a modulo q map 

Mod, : ^ (M/gZ)^, (Mod<,/)(x) = f{x) (mod g), / e R^. 

Then Mod,(Z^) = (Z/qZ)^ is a subgroup of the toric group (R/gZ)^, and MoAq{F{G, e; Z)) 
is a subgroup of Mod5(Z^). 
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Given orientations p,(T & 0{G); there is an involution Qp^a- ■ [0,q]^ [0, 9] 
defined by 



where g G [O,*;]^ and v is an end-vertex of the edge e. Clearly, Q^^^ is a bijection 
from [0, q]^ to [0, q]^, and is also a bijection from (0, q)^ to (0, q)^ , where (0, q) = 

{X eR\0 < X < q}. Moreover, Qp^aQa.e = Qp.E- 

Recall that two orientations £1,62 G 0{G) are said to be Eulerian equivalent, 
written £1 ~ £2, if the induced spanning subdigraph by the edge subset 

£'(£1 7^ £2) ei{v, e) ^ S2{v, e),v is an end- vertex of e} 

is a directed Eulerian subgraph with the orientation either £1 or £2. 

Lemma 5.1. (a) The relation ~ is an equivalence relation on 0{G). 

(b) Let p, cr S 0{G) be Eulerian equivalent. If p is totally cyclic, so is a. 

(c) Let p,a € 0{G) and p ^ a. Then Qp^a '■ Q^thi^^^^) ~^ Q^thi^^p) ^ 
bijection, sending lattice points to lattice points. In particular, 

Qp,.{qA+{G,a)) = qA+{G,p), 

iPp{G,q) = ipa{G,q), 

Cpp{G,e;q) = Cpa{G,q). 

Proof, (a) The rcflexivity and the symmetric property arc obvious. Transi- 
tivity is a straightforward computation. Note that a digraph [H, p) is Eulerian if 
and only if for all v e V{H), 




Let £i e 0{G) {i 



1, 2, 3) be such that £1 ~ £2 and £2 ~ £3. Then 



^ £2(u,e)= ^ £2(w,e) = 0. 



Thus 



E 



si{v,e) 



^ ei{v,e)+ ^ ei(i',e) 



eeE{ei^e3) 



eeB(£i=£25^e3) e6-E(ei#e2=e3) 



eeE(£i=e2#e3) eeE{sijte2=S3) 



E S2{v,e) 



eeE(£i=£25^£3)U-E(ei#e2#£3) 



ee£;(£i5^£2=£3)U£;(£l5<i£2#£3) 



Y £2(v,e)- Y £2(^,6) = 0. 



eeE(£27^£3) eeB(£i5^£2) 



This means that £1 is Eulerian equivalent to £3. 



14 



BEIFANG CHEN AND RICHARD STANLEY 



(b) Suppose {G,p) contains a directed cut {U,eu), where U = [S,S'^] with 
5 C y. Since E{p ^ a) is directed Eulerian with the orientation p, then by 
Proposition I3.1f b) , 

^ [p,eu]{x)=0. 

Since [p, eu] = 1, it fohows that U r\E{p / cr) = 0. This means that e) — <j{v, e) 
for all edges e e C/, where v is an end- vertex of e and v £ S. So {U, Su) is a directed 
cut of (G, e). This is a contradiction. 

(c) For a flow / G (7A+ (G, a) (/ £ (zA+ (G, cr)), we have 

X! P{'"^e.)i.Qp,<yf){(^) = X! '^(^'e)/(e) +g ^ p(w,e) = 0. 

This shows that Qp^cr/ £ 9A+(G,p) {Qp^f G 9A+(G, p)). Clearly, Qp^o- sends 
lattice points to lattice points by definition. Therefore, (pp{G^q) = (pa{G^q) and 
'/fp(G, g) = (pa(G, g). □ 

For two Eulerian equivalent orientations p, cr, we have seen that the digraph 
(G, p) contains no directed cut if and only if (G, cr) contains no directed cut. This 
means that ^ induces an equivalence relation on Otc(G); and each equivalence 
class of ~ in C'tc(G) is an equivalence class of ~ in 0{G). We denote by [C'tc(G)] 
the quotient set ©^(G)/^ of Eulerian equivalence classes. For each p G CtcCG), 
let [p] G [0Tc(G)] denote the equivalence class of p, and define 

<P[p] (G, q) = <^p(G, q) and (^[pj (G, g) = (^p(G, g). 

The following nontrivial Lemma is due to Tutte. It is crucial to the proof of our 
main result Theorem II. 2 1 so we present a proof with our notations. 

Lemma 5.2 (Tutte [23]). The map Mod, : F^{G,e;q) -> F{G,e;Z/qZ) and its 
restriction Modg : Fnz^{G,£;q) — > i^nz(G, e; Z/qZ) are surjective. 

Proof. The second part of the lemma implies the first part. In fact, every 
flow / G F{G, e; ZgZ) can be viewed as a nowhere-zero flow f\Ef on the subdigraph 
{V,Ef,e), where Ef := {e G E \ /(e) ^ 0}. Let g be a nowhere-zero integer g-flow 
on {V,Ef,e) such that Modq{g) — f\Ef Then g is extended to an integer g-flow 
on (G, e) by setting g = l on E — Ef. 

Now for each / G F^{G, e; q) we write / = Modg/. We identify Z/gZ with the 
set {0, 1, . . . , g — 1} and view each modular flow g G F{G, e; Z/gZ) as an integer- 
valued function g : E {0, 1, . . . , g — 1}. Then the map Qp^^ maps F{G, e; Z/gZ) 
to F(G,p;Z/gZ), and Fnz^(G, e; Z/gZ) to F^^^iG, p;1/ q1). For each g G [Q,q]^ 
and p G 0{G), we define 

Fix a modular fiow / G i^nzz:(G, e; Z/gZ). Let p* be a particular orientation on G 
such that 

v{Qp'J:P*) - min{rKQp,e/,p) | P G 0(G)}. 
We write /* Qp*,,./ and define / := Pe,p' f* ■ 
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If r;(/*,p*) = 0, then EeeB/5*(«'e)/*(e) = for all v e F. This means that /* 
is an integer q-flow of (G, p*). Whence / is an integer g-flow of (G, e). By definition 
of Pe,p* and Qp',e, we see that 



/(e) = 



/(e) if e(w,e) = p*(u,e) 
f{e)-q a e{v,e) p*{v,e) 



= /(e) (mod q) for all e e E. 

The surjectivity of Modg : Fnzz{G,e;q) Fnz{G,e;Z/qZ) follows immediately. 

We now claim that r]{f*,p*) — 0. Suppose ri{f*,p*) > 0. Then there exist a 
vertex u d V and an integer k such that 

J2p*i^,e)r{e) = kq>0 {<0). (5.1) 

eeE 

Note that for any function g : E ^ A, where A is an abelian group, and for any 
orientation p G 0{G), we have 

^^p(t;,e)5(e) = 0. 

vev eeE 

In particular, for the function /* and the orientation p* , we have 

^^p*(^;,e)/*(e)=0. 

Since ()5.ip and /* is a modular fiow, there exists a vertex w such that 
Y,P*iw,e)rie)^-lq<0 (>0). 

eeE 

Thus there is a path P — vqCiVi ■ ■ ■ CnVn with vq — u and w„ = w, such that 
P*{vo,ei) = 1, p*(w„,e„) = -1, and/9*(?;i,ei)p*(?;i,ei+i) = -1, where 1 < i < n-1. 
Let p** be an orientation on G given by 

-p*{v, e) if e = Ci for some 1 < « < ri, 
p*(u,e) otherwise. 



P (w,e) 



We write /** := Qp",^./*. Then /** = Qp",p.Qp.,./ = Qp",J. 
Notice that at each vertex v ^ V, we have 

5]p"(«,e)/**(e) = ^p*(t;,e)/*(e)-g ^ p*iv,e). 

eeE eeE eeE{p*^p") 

In particular, for the vertices u,w, and other vertices v different from u and w, we 
have 

^ p**(^., e)/**(e) = ^ p*(^., e)/*(e) - g - (fc - l)g, 

eeE eeE 



^ p**(u;, e)/**(e) = ^ p*(«;, e)/*(e) + g = (1 - 0?, 

eeE eeE 

^p**(z;,e)/"(e) = ^p*(«,e)/*(e). 



eeE eeE 

It follows that 

'7(/**,P**)='7(/*,P*)-2<z <??(/*, p*). 
This is contradictory to the selection of p* that ??(/*, p*) is minimum. □ 
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For each real-valued function / : _B — s- M and any orientation p on G, we 
associate with / and p an orientation p/, defined for each {v,e) x E by 

For two orientations p, cr G 0{G), we associate a symmetric difference function 
Ip^cr : i? — )■ {0, 1}, defined for each edge e e i? (and its one end-vertex v) by 

/p,.(e) = ( ? ^[^[^'^jT'^J^'^j' (5.3) 
[1 it p('y, e) (t(w, e). 

Lemma 5.3. Let /i,/2 £ F{G,e;q). If fi{e) = /2(e) (mod g) for all edges e E E, 
then e and e are Eulerian equivalent. 

Proof. Let us simply write :=£/., i = 1, 2. It suffices to show that at each 
vertex i), 

£2(«,e)=0. (5.4) 

eG-E(ei5^e2) 

Since /i(e) = /2(e) (mod g) for all e £ E, then 

/i(e) = /2(e) + fleQ, e e 
where G { — 1,0,1}. More precisely, for each edge e at its one end- vertex v, 

{/2(e) if ei(w,e) = £2(w,e), 

/2(e) - g if ei{v, e) 7^ £2(w, e) = e{v, e), 
/2(e) +q if ei(w,e) 7^ £2(w,e) 7^ e{v,e). 

Let = Pg.^i^fi, i — 1,2. Then gi are real g-flows of {G,ei), and for an edge e at 
its one end- vertex v, 

, , f g2(e) if £i(w,e) = £2(-y,e), 

' 1 9-32(e) if £i(w,e) 7^ £2(w,e). 

Note that at each vertex v £ V, 

^£,(z;,e)<7.(e)=0, i = l,2. (5.5) 

eG-B 

Now consider the case of i = 1 in (|5.5p . Replace 51(e) by 52(e) if £i(u, e) = £2(w, e) 
and by q — 52(e) if £i(w, e) 7^ £2(w, e); we obtain 

£2(w,e)52(e) - ^ £2(-u,e)(g- 52(e)) 
= ^£2(w,e)52(e) ^ £2(u,e)=0. 

eG-E ee_E(£i#e2) 

Applying (|5.5p for i 2, we see that (|5.4p is true. □ 

Lemma 5.4. Let p, cr, a; G Otc(G) &e Eulerian equivalent orientations, and let 
f G qA'^i^{G, e) be a real q-flow. Then 

(a) Sf = p. 

(b) Pe,.Q.,p-Pp,e(gA^;aG,£)) =(ZA?,(G,£). 

(c) Pe,aQcr,pPp,ef = Pe,u,Qi^,pPp,ef «/ a?^C^ 071^?/ if a = UJ . 

(d) F(G,£;g) nMod-l(Mod,/) = {Pe,aQc.,pPp,ef \ a ^ p} . 
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Proof, (a) By definition of A^^{G,e) and the fact / G qA^i^{G,e), we have 
[p, e](e)/(e) > for all e & E. So for each edge e at its one end- vertex v, p{v, e) — 
e{v,e) if /(e) > 0, and p{v,e) ~ —e{v,e) if /(e) < 0. By definition of Sf, we see 
that £f = p. 

(b) Recall Lemma 14.21 and Lemma IS.lT c) that 

Pp,,{qAP,{G,e)) = qA+,{G,p) and Q„,p{qA+{G,p))=qA+,{G,<j). 

Since Pg^a is an involution, then Fe^o^(gA+ (G, ct)) = qA^^{G,e). The identity fol- 
lows immediately by composition. 

(c) The sufficiency is trivial. For necessity, we write g :~ Pe,aQa,pPp.e.f and 
h ■= Pe,ujQuj,pPp.Ef ■ Since g S qA'^^{G,e) and h £ qA'^^{G,e), we have Sg = cr and 
Eh — u!. Clearly, ii g — h, then a — w. 

(d) Let us write g := Ps,aQa,pPp,ef for SiH orientation a such that a ~ p. By 
definition of Pe,a,Qa.,p, Pp,£, we have g G A"^{G, e) and 

^(e) = /(e) -f flgg, eeE, 

where S { — 1,0,1}. Clearly, g{e) = f{e){inod q) for all e £ E; namely, g G 
Mod^^(Modq/). Conversely, let h G F{G,e;q) be such that Modqh = Mod,/. We 
have 

h{e) = f{e) + beq, eeE, 
where be G { — 1,0, 1}. By definition of e/,^^, a straightforward calculation shows 
that for each non-loop edge e (and its one end- vertex v), 

/(e) if h{e) > 0, /(e) > (^ £/i(w, e) = £/(i;, e) = £(t;, e)), 

/(e) if /i(e) < 0, /(e) < (^ s,,{v, e) = Sfiv, e) ^ £(t;, e)), 

/(e) - g if h{e) < 0, /(e) > (^ £,,(«, e) =^ £/(t;, e) = £(w, e)), 

/(e) + g if h{e) > 0, /(e) < (^ s,,{v, e) ^ £/(i;, e) ^ £(t;, e)). 

By definition of P^sh ? Qch^sf i another straightforward calculation shows that 

for each non-loop edge e (and its one end- vertex v), 

if £/i(w, e) = £/(w, e) = £(w, e), 
if £/i(w, e) = £/(w, e) 7^ £(w, e), 
if eh{v, e) 7^ £/(i;, e) = e{v, e), 
if £h(i;, e) 7^ £/(i;, e) 7^ e{v, e). 

This means that h = P^^g^f^Q^^^^jPg^j^^f. Since £/ = p by Part (a), thus h = 

Pe,aQa,pPp,sf with Q; = £ft. □ 

Proposition 5.5. T/ie number of orientations on G that are Eulerian equivalent 
to £ is the number of O-l flows of (G,e), i.e., 

#[£] - (^e(G, 1) = |A+ (G,£) n Z^|. (5.6) 

Proof. It is enough to show that the following map 

[£] = {p G 0{G) \p^e}^ A+ (G, £) n Z^, p V, 

is a bijection. Fix an orientation p that is Eulerian equivalent to e. Note that Jp^^ is 
the characteristic function of the edge subset E[p 7^ £). Since E(p 7^ £) is directed 
Eulerian with the orientation e. Then /p^^ is a flow of (G,£) with values in {0, 1}. 
So the map p f-> /p.g is well-deflned, and is clearly injective. Conversely, given a 0-1 
flow / of (G, e). Let p be an orientation on G defined by p(i;, e) = £(z), e) if /(e) = 
and p{v,e) — ~£{v,e) if /(e) — 1, where v is an end-vertex of the edge e. Then 



h{e) 



(P,,e,Qe,^e,Pe,,./)(e) = 
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E{p ^ e) = {e G E{G) \ /(e) = 1}, which is directed Eulerian with the orientation 
e. This means that p ^ e and Ip^e = f ■ Hence the map is surjective. □ 

Proof of Theorem 11.21 

It has been shown that (p{G,q) is a polynomial function of degree n{G). Fix 
an orientation p G ©^(G). For each orientation a ^ p and any / G qA^^{G, e), we 
have 

|F(G,e;<7)nM0d-l(M0dJ)| = #{Pe,aQa,pPp,e/ I « ^ ^} 

- #M=#[P]. (5.7) 
Now apply Parts (b) and (d) of Lemma 15.41 we have the disjoint unions 

cre[p] crelp] 

□ {Pe,.Q.,pPp,ef} 
<T6[p],/e?Af,(G,e) 

= □ F(G,£;g)nMod-\Mod,/) 

= F{G,e;q)nMod-^Mody{qAP^{G,e)). 

Note that the above orientation p can be replaced by any orientation a that is 
Eulerian equivalent to p. We further have 



□ gA-,(G,£)=F(G,e;g)nMod-iModJ □ <zA^,(G,£; 
ae[p] ^c^e[p] 



(5.8) 



On the one hand, recall Lpa{G, q) — ipp{G, q) whenever a ^ p (see Lemma lS.ip : then 
the number of lattice points in the left-hand side of (|5.8p is 

Vp{G,q)-m- 

On the other hand, (15. 7p implies that the number of lattice points of the right-hand 
side of is 

Mod,( □ gA-,(G,£)nZ^)| 

cre[p] 

It then follows that 

Mod,( □ (7A^,(G,e)nZ^)| =^p(G,g). 

o-G[p] 

Note that (14. 3|) implies the disjoint decomposition 

F,.UG,e-q)^ U □ gA,-,(G,e)nZ^. 

pe[0^c(G)\ a£[p\ 

Applying the map Modg, we obtain the following disjoint decomposition 

Fn,(G,£;Z/(7Z) = Mod, (Fn,,(G, £;(?)) 

- □ Mod J □ gA^,(G,£)nZ^ 

[ple[OTc(G)] ^<TG[p] 



(5.9) 
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The first equation follows from the surjectivity of Modq by Lemma 15.21 To see 
the disjointness of the miion in the second equation, suppose the union is not 
disjoint. This means that there exist integer flows fi G qA'^l{G, e) and orientations 
Pi such that (Ti Pi {i = 1,2), [pi] ^ [p2], and Modg(/i) = Modg(/2). Then 
Lemma I5.4r a) implies that Sf- = ai, and Lemma 15.31 implies that £/i £/2- It 
follows by transitivity that pi ~ p2, i-e., [pi] = [p2]- This is a contradiction. 
Counting the number of elements of both sides of (|5.9p . we obtain 

^{G,q)= fpiG,q)- 
[p]e[OTc(G)] 

The Reciprocity Law ()1.12[) follows from the Reciprocity Law (|1.7p and the defini- 
tion of <f>{G, q). 

Let q^O. We have ip{G,0) = (-1)"(^)(^(G, 0) by (frT2]) . Since ^p(G,0) = 1 
for all p G Otc(G), we see that (^(G,0) = #[Otc{G)] by (lllill . Hence ip{G,0) = 
{-iy(G)^[0,,XG)]. □ 

Searching online we found a result on modular flow reciprocity by Breuer and 
Sanyal (2, which states that (— 1)"('^V(G, — g) counts the numbers of pairs (/, cr), 
where / is a flow of (G, e) modulo q and is a totally cyclic reorientation of the 
digraph (G/supp/,e). The result can be written as the sum 

V{G,-q) ^ (-1)"('^) ^m,q)\0,c{G/X)l (5.10) 

XCE 

where G/X is the graph obtained from G by contacting the edges of X. The term 
ip{{X),q) \Otc{G/X)\ in (fO))) is nonzero if and only if the graphs (X) and G/X 
are bridgeless. The formula (|5.10p can be argued straightforward as follows. 

Notice the trivial fact that each flow / corresponds to a nowhere-zero flow on 
its support supp (/) := {e G I /(e) 7^ 0}. This means that 

YCX 

The Mobius inversion implies 

^{{x),t)^ ^(-i)i^-^ir<^>, XCE. 

YCX 

In particular, for X ~ E and t = —q, we have 
^(G,-g) = 

Y<ZE 

- 5:(-i)i^-^i+"<^)5:^((x),,) 

Y<ZE XCY 

= E (-i)i^-^i+"<^>. 

XCE XCYCE 
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Since n{Y) ^ n{X) + n{{Y) / X) for all edge subsets F of G such that X C y, we 
see that for each fixed edge subset X C E, 



E = E (-1)' 



^^\E-Y\+n{X)+n({Y)/X) 

XCYCE XCYCE 

^ (^_l-^n{X) ,^_^-^\E{G/X)-Z\+n{Z) 



ZCE{G/X) 

{-l)^(x)+^iG/x)\o^^,(^G/X)\. 



The above last equality follows from the Zaslavsky formula p.4p about the flow 
arrangement ^FL(G,e;M), i.e., 

^(G, -1) = E = (-ir(«)|o,c(G)|. 

ZCE 

Now the identity (IS.lOp follows immediately. 

6. Connection with the Tutte polynomial 

The Tutte polynomial (see [3], p. 337) of a graph G = (V, E) is a polynomial in 
two variables, which may be defined as 

TG{x,y)^ J2{^~lY^''^-'''^^Hy-ir'^^\ (6.1) 

ACE 

where (E) = \V\ - c(G), r{A) = \V\ - c{A), n(A) ^ \A\ - r{A}. The polynomial 
Tcix^y) satisfies the Deletion-Contraction Relation: 

{xTQ/e[x, y) if e is a bridge, 

yTG-e{x, y) if e is a loop, 

Tc-e {x,y) + Ta/e {x,y) otherwise. 

It is well-known that the flow polynomial ip{G,t) is related to TG{x,y) by 

^(G,t) = (-l)"(«)TG(0,l-i). (6.2) 

Thus 

^(G, t) = (-!)"(«) (p(G, -t) = Tg(0, t + 1). 
We conclude the information as the following proposition. 

Proposition 6.1. The Tutte polynomial Tq is related to ip and (p as follows: 

Tg{0, t) = (^(G, t-\)^ (-1)"(«V(G, l~t). (6.3) 

In particular, Tg{0, 1) = |(/3(G,0)| = (^(G, 0) counts the number of Eulerian equiv- 
alence classes of totally cyclic orientations on G. 

Example 6.2. Let Bn(u, v) be a graph with two vertices u, v, and n multiple edges 
ei, . . . , e„ between u and v. Let e be an orientation of Bn(u, v) such that all edges 
have the direction from u to v. The number of integer flows / of Bn{u, v) such that 
I/I < g is equal to the number of integer solutions of the linear inequality system 

xi ^ \-Xn = 0, -{q - I) < Xi < q - 1. 

Let Xi — yi — (q — I) . The above system reduces to 

yi + --- + y„ = n(g-l), 0<y. <2g-2. (6.4) 
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Recall that the number of nonnegative integer solutions of yi + • • • + j/n = r is 

'n + r^l\ /ri + r — 1' 
/ \ n — 1 

So the number of nonnegative integer solutions of j/i + ■ • ■ + y„ = n{q — 1) is 
( ^n-i ) ' "^-^lich includes the number of integer solutions of (|6.4I) and the number of 
nonnegative integer solutions having at least one Ui > 2q — 1. 

Let Sn{q) denote the number of integer solutions of (|6.4|) . To figure out Sn{q), 
let Y be the set of nonnegative integer solutions oi yi + ■ ■ ■ + yn — n{q — 1), 
Yi the set of nonnegative integer solutions of j/i + • ■ • + y,i = n{q — 1) with the 
ith variable yi > 2q — 1, and Yq the set of integer solutions of (|6.4I) . Then Yq = 
Y — [j"^^ Yi. Consider the case where j variables are greater than or equal to 2q — l, 

say, yn-j+i,- ■ ■ , yn] then yi H \-yn-j < n{q -1) - j{2q-l). The number of such 

solutions is equal to the number of nonnegative integer solutions of the equation 
2/1 H 1- 2/n-i + y' = n{q - 1) - j{2q - 1), which is given by 

j + 1 \ / {n-2j)q \^/{n-2j)q' 

n{q-l)-j{2q-l)/ {n{q - 1) ~ j{2q ^ 1) J { n-j 
Applying the Inclusion-Exclusion Principle, 

07^7C[n] ^iel 

We list si,S2,S3,S4 explicitly as follows: 

si(g) = 1, S2{q)=2q-l, 53(5) = Sg^ - + 1, 

sM = (^'^3" - 4 (^^g'^) = i(2<z - 1) (8g^ - 89 + 3) . 

Let So{q) (= 1) denote the number of zero flows of Bn{u, v). Then ip^{Bi, q), defined 
as the number of nowhere-zero integer flows of Bi{u,v), counts the number of 
integer flows / such that < \ f\ < q, and is given by 

^.{B4,q) = #([l-q,q-l]''^'''^nF{B4,e;Z)- [j F^ (6.5) 

^ e(iE{Bi) ' 

= 54(9) - 4s3(g) + 6s2(g) - 4si(g) + so(g) 
= ^(g-l)(8g2-22g + 2l). 

Likewise, the number of flows of Bn{u,v) modulo q is <?"^^, n > 1. Thus 
ip{B4.,q), defined as the number of nowhere-zero flows of Bi{u,v) modulo q, is 
given by 



v{B4,q) = if(F{B4,e;Z/qZ)- \J Fe] 



(6.6) 



eeE(B4 

q^ - 4q^ +6q- 4 + 1 
(g-l)(g2-3g + 3). 
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In particular, 

1^,(^4, 0)1 = \^{Bi,-l)\ = 14, |(^(B4,0)| =3. 

There are 14 totally cyclic orientations on B4{u,v) by Theorem 11.11 The 14 ori- 
entations can be grouped into 3 Eulerian equivalence classes by Theorem II. 2t see 
Figure [H 




Figure 1. The three Eulerian equivalence classes of Otc{B4). 

Let Sk be an orientation on Bn(u, v) such that the arrows of the edges ei, . . . , 
point from u to v, and the arrows of the edges Ck+i, . . . , e„ point from v to u. If a 
0-1 flow / of {Bn,£k) has value 1 on exact j edges of ei, . . . , e^, then / must have 
value 1 on exact j edges of e^+i, ■ ■ ■ , e„. Thus by Lemma 15. 5[ ip^k {Bn, 1) is the 
number of 0-1 flows of {Bn,Sk), and is given by 

min(k,n — k) 

E ( )("-')■ 

For the case n = 4, we see that (peo[Bi,l) = 1, Lp^^[Bi,l) = 4, (^^^{Bi,!) = 6, 
(fe^iBi, 1) — 4, and (pe^iB^, 1) — 1. According to Theorem 11.21 "^£^(^4, 1) counts 
the number of orientations of B^^ that are Eulerian equivalent to Sk ■ Indeed, there 
are 4 orientations Eulerian equivalent to £i and £3 respectively, and 6 orientations 
Eulerian equivalent to £2- These orientations are listed in Figure [TJ However, the 
Eulerian equivalence classes for Eq and £4 are singletons; the integral and modular 
flow polynomials with respect to these orientations are the zero polynomial. We 
list the two orientations as follows: 




Figure 2 . The other two Eulerian equivalence classes of orienta- 
tions of S4 with directed cut. 
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Remark. The coefRcients of the integral flow polynomial (G, t) are not neces- 
sarily integers as shown in Example 16.21 The combinatorial interpretation on the 
coefRcients of ip^{G,t) is particularly wanted. 
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